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Abstract 


The  problem  considered  here  is  that  of  a  finite,  riaid  insert  partially 
embedded  in  and  adhesively  bonded  to  an  elastic  naif  plane.  Two  distinct 
problems  are  investigated:  the  shear  pullout  of  the  insert  without  rotation, 
which  takes  into  account  the  adhesive's  resistance  to  shear  deformation,  and 
an  opening  problem  which  incorporates  the  adhesive's  resistance  to  normal 
deformation.  This  latter  problem  assumes  the  presence  of  an  edge  crack  in  the 
half  plane  subjected  to  opening  pressure  which  equals  in  magnitude  the  normal 
stress  distribution  due  to  the  pullout  problem.  These  mixed  boundary  value 
problems  are  governed  by  singular  integral  or  integrodi fferential  equations. 
Numerical  results  are  obtained  via  a  technique  Introduced  by  Gerasoulis  and 
Srivastav.  Several  important  physical  quantities  are  calculated,  such  as  the 
shear  and  normal  stress  distributions  along  the  bonded  interface,  and  the 
crack  opening  displacements. 


1.  Introduction 


The  problem  of  a  rigid  insert  partially  embedded  in  and  adhesively  bonded 
to  an  elastic  half  space  is  considered  within  the  context  of  plane  elasto- 
statlcs.  The  thickness  of  the  insert  Is  assumed  to  be  negligible.  It  is  also 
assumed  that  the  adhesive's  deformation  remains  in  the  linear  elastic  range. 
The  thickness  of  the  adhesive  is  small,  but  not  negligible. 

Loading  is  applied  to  this  insert  so  that  it  translates  in  the  vertical 
direction  without  rotation  (Fig.  1).  Of  interest  here  are  the  shear  and 
normal  stresses  generated  along  the  bonded  interface.  These  stresses  will 
cause  the  adhesive  to  undergo  both  shear  and  normal  deformation.  Two  distinct 
problems  are  considered:  the  shear  pullout  problem  already  described,  which 
takes  Into  account  the  adhesive's  resistance  to  shear  deformation,  and  the 
opening  problem,  described  next,  which  considers  the  adhesive's  resistance  to 
normal  deformation.  The  opening  problem  assumes  an  edge  crack  in  a  half  space 
subjected  to  opening  pressure  which  1$  equal  In  magnitude  to  the  normal  stress 
distribution  resulting  from  the  pullout  problem.  In  this  case,  the  adhesive 
material's  resistance  to  normal  deformation  Is  used  to  reduce  the  opening 
pressure's  Intensity.  The  physical  quantities  of  Interest  in  this  problem  are 
the  crack  opening  displacements. 

The  problems  considered  here  can  be  viewed  as  appropriate  idealizations 
for  studying  two  separate,  yet  similar,  classes  of  problems.  In  the  context 
In  which  they  have  been  presented  thus  far,  they  are  problems  in  plane 
strain.  As  such,  they  are  appropriate  for  investigating  the  mechanics  of 
Interface  failure  In  fiber-reinforced  composite  materials.  The  more  realistic 
case  of  adhesive  failure  (adhesive  material  allowed  to  deform  non-linearly)  is 
currently  being  considered  by  the  authors.  As  pointed  out  by  Brussat  and 


Westmann  [1]  and  by  Chamis  [2J,  Interface  damage  is  believed  to  be  one  of  tne 
earliest  forms  of  damage  in  such  composite  materials.  The  subject  of  fiber 
debondiny  has  attracted  several  investigators  over  the  last  several  years  (see 
Cl],  [2]  for  summary).  One  fundamental  question  that  seems  to  remain  still 
unanswered,  and  which  is  addressed  in  this  paper,  is  whether  tne  surfaces  of  a 
crack  at  the  fiber-matrix  Interface  are  traction- free,  or  remain  in  contact. 

It  is  suggested  in  [U  that  the  debonded  surfaces  may  be  open  or  closed 
depending  on  the  relative  magnitude  of  the  applied  versus  the  residual 
stresses. 

In  the  generalized  plane  stress  context,  these  problems  deal  with  finite 
rigid  inserts  of  negligible  width,  partially  embedded  within  a  semi-infinite 
elastic  sheet.  In  this  case,  the  results  obtained  may  serve  to  better  under¬ 
stand  the  stress  distribution  at  the  bonded  Interface  of  dissimilar  materials, 
where  one  is  significantly  more  rigid  than  the  other.  More  importantly,  an 
extension  of  this  work  which  will  be  addressed  in  the  near  future  will 
consider  the  case  of  an  elastic  insert  of  finite  width  embedded  in  an  elastic 
half  plane.  That  problem  would  be  applicable  to  several  physical  situations, 
such  as  the  primary  adhesively  bonded  aircraft  fuselage  being  considered  by 
the  U.S.  Air  Force  [3]. 

The  problems  are  formulated  as  mixed  boundary  value  problems.  This  leads 
to  singular  Integral  or  Integrodl fferentlal  equations,  which  are  solved  numer¬ 
ically  for  the  desired  physical  quantities  by  use  of  a  numerical  technique 
Introduced  by  Gerasoulls  and  Srlvastav  [4],  [5].  The  material  constants 
selected  for  the  adhesive  and  for  the  adherend  material  correspond  to 
materials  commercially  available. 
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2.  Formulation  and  Basic  Equations 

The  formulation  of  the  problems  considered  in  this  investigation  makes 
use  of  certain  results  obtained  in  two  earlier  papers  [6],  [7].  Those  papers 
considered  the  problem  of  a  finite,  rigid  block  embedded  in  an  elastic  half 
space  [6],  and  that  of  a  rectangular  trench  in  an  elastic  half  plane  [7].  The 
problems  were  formulated  by  superposition  of  the  solutions  for  rigid  line 
inclusions  and  cracks,  respectively,  parallel  and  perpendicular  to  the  free 
surface  of  the  half  plane.  The  results  which  pertain  to  the  present  work  are 
those  derived  in  [6]  for  a  rigid  line  inclusion  perpendicular  to  the  free 
surface  of  the  half  plane,  and  the  ones  derived  in  [7]  for  an  edge  crack. 

A.  Rigid  Insert  Pullout:  The  Shear  Node 

A  rigid  line  inclusion  is  partially  embedded  in  and  bonded  to  an  elastic 
half  plane,  y>0,  so  that  it  occupies  the  line  segment  0<y<h,  x=0  in  the  half 
plane.  The  Inclusion  Is  loaded  by  a  force  P  acting  along  the  negative  y-axls 
as  shown  in  Fig.  1.  The  material  properties  of  the  half  plane  are  taken 
as  w  and  u  is  the  shear  modulus,  while  k  Is  related  to  Poisson's 
ratio  v  by  k  ■  3-4v  (plane  strain),  or  k  »  (3-v)/(l+v)  (generalized  plane 
stress).  The  Inclusion  Is  assumed  to  be  bonded  to  the  half  plane  by  use  of 
some  adhesive  material  In  such  a  manner  that  there  are  no  discontinuities  In 
the  application  of  the  adhesive  to  the  Inserted  inclusion,  and  no  significant 
variations  In  Its  thickness,  which  is  also  assumed  to  be  negligibly  small 
(less  than  5x10"^  in.).  The  adhesive  Is  assumed  to  behave  as  a  number  of 


linear  springs,  subjected  to  shear  in  this  case.  The  stiffness  of  these 
springs,  ks,  is  computed  using  the  mechanical  properties  of  the  adhesive. 


Tnis  is  discussed  in  detail  in  a  later  section. 

The  shear  stress  developed  in  the  adhesive  (and  transferred  to  the  half 
plane)  may  be  written  as  follows: 


Txy(x,y)  =*  ks  [u  (x,y)  -  uQ]  ;  x  *  0,  0  <  y  <  h 


(2.1) 


Here,  u0  is  a  constant  representing  the  y-di spl acement  of  the  rigid  inclusion, 
while  uy  designates  the  y-displacement  of  the  half  plane  material.  Thus,  the 
quantity  in  the  brackets  gives  the  displacement  of  the  insert  relative  to  the 
half  plane,  which  is  also  the  shear  deformation  of  the  adhesive.  The  constant 
ks  is,  as  defined  earlier,  the  stiffness  of  the  adhesive  modeled  as  a  shear 
spring. 

The  displacement  uQ  cannot  be  calculated.  Differentiation  of  (2.1)  with 
respect  to  y  yields  the  following  equation: 


lyOgU.y)]  •■^|?Ctxy(x.,)3  ;  x  -  o. 


(2.2) 


The  right-hand  side  of  eqn  (2.2)  is  given  in  terms  of  a  singular  integral 
equation  derived  in  [6]  using  Integral  transform  techniques  [8],  [9].  When 
that  result  Is  used,  equation  (2.2)  takes  on  the  following  form: 


4w(  K+l  ) 


h  P(t)dt 
o  y-t 


dt  h  nr \  r  Lsili  2[2t-ic(y+t)  ] 
t  o  L  y+t  ,_*,2 


(y+t)‘ 


(y-t)' 


]  dt)  *  -r—  Cr  (x,y)]  ;  x  =■  0,  0  <  y  <  h 


*s  3y  xy 


(2.3) 


The  function  Q(y)  is  the  unknown  shear  stress  discontinuity  associated  with 
the  rigid  inclusion  and  it  is  defined  as  follows: 
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0  ( y )  3  T 


(2.4) 


xy 


x>0,  0<y<h 


The  superscripts  (1)  and  (2)  denote  the  regions  to  the  left  and  to  the  riyht 

of  the  inclusion,  respectively.  The  symmetry  of  the  problem  requires 

that  t. .  ^  3  -t  Clearly  then,  the  following  relate  must  hold 

-y  xy 

between  T  and  D: 

*y 


txy(0,y)  3  D(y)  ;  0  £  y  _<  h 


(2.5) 


Substitution  of  eqn  (2.5)  into  eqn  (2.3)  leads  to  the  following  integrodif- 
ferential  equation  in  the  unknown  D(y) : 

/J  +  /J  0(t)K(y,t)dt  3  -  |y  CD(y)]  ;  x=0,  0  <  y  <  h  (2.6) 

Here,  the  function  K(y,t)  is  given  as  follows: 

K(y.t)  3  -  i-  [  3  (2.7) 

Zk  (y+t)Z  (yt)3 

It  should  be  noted  here  that  the  expression  for  |-  u  substituted  in  eqn  (2.2) 

ay  y 

to  obtain  eqn  (2.3)  was  derived  in  [6]  so  that  the  following  conditions  are 
satisfied: 

„<*>  .  „<»)  .  o 

X  X 

}  x-o,  0  <y  <  h  ;  (2.8) 


Tyy  =  Txy  =  0  *  y  =  °*  0  £  lxl  < 


Equation  (2.8)  requires  that  the  displacements  on  either  side  of  the  inclusion 
be  continuous,  while  eqn  (2.9)  clears  the  half  plane  surface  of  stresses. 

Equilibrium  in  the  y-direction  is  satisfied  by  requiring  that  the  follow¬ 
ing  relationship  hold: 


/"  0(y)  dy  =  P 


(2.10) 


Thus,  the  solution  of  the  shear  pullout  of  the  rigid  insert  problem  consists 
of  solving  the  governing  integrodi fferential  equation  (2.6)  for  the  unknown 
D(y)  subject  to  the  equilibrium  condition  (2.10).  The  solution  is  obtained  by 
means  of  numerical  methods  as  will  be  discussed  in  the  numerical  analysis 


section. 


Once  D(y)  has  been  determined  the  stresses  generated  along  the  bond  line 
are  readily  obtained.  The  shear  stresses  are  computed  usiny  eqn  (2.5),  while 
the  normal  stresses  are  obtained  from  the  following  [6]: 


Txx{0’y)  *  17&YT  {/5  +  Jo  D(t)  L(y’t)  dt}  :  0  y  -  h  (2*u) 


The  function  L(y,t)  Is  defined  as: 


1  r  3c*l  .  2( 3t-iey)  8yt  i 


B.  Rigid  Insert  Pullout:  The  Opening  Node 


The  normal  stress  distribution  generated  along  the  bond  line  of  the 
insert  during  pullout  is  now  applied  as  opening  pressure  on  an  edge  crack  in 
the  half  plane.  It  Is  assumed  that  the  crack  extends  from  y= 0  to  y=h  and  it 
is  located  at  the  bond  line,  x=0.  The  adhesive  material  in  this  case  behaves 
as  a  number  of  linear  springs  of  stiffness  k  subjected  to  tension.  Thus,  it 
resists  the  opening  of  the  crack  by  providing  a  stress  equal  to  k(C0D/2)  which 
reduces  the  intensity  of  the  applied  opening  pressure.  The  crack  opening 
displacement  (COO)  is  the  distance  between  the  crack  surfaces. 

The  normal  stresses  within  a  half  plane  containing  an  edge  crack  are 
given  in  [7].  The  bond  line  stresses  are  given  in  terms  of  the  dislocation 
density  B(y)  as: 


Txx(0,y) 


h  B(t)dt 
o  t-y 


+  /S  B(t)  [ 


Wf'p  ]  dt}  ; 
(yn)3 


0  <  y  <  h 


(2.13) 

Here,  u  and  <  are  the  half  space  material  constants  given  earlier.  The  super¬ 
script  c  is  used  to  distinguish  these  normal  stresses  from  the  ones  associated 
with  the  shear  pullout  problem.  The  dislocation  density  B(y)  is  defined  as: 


B(y)  [u(*>  -  ;  x  -  0,  0  <  y  <  h  (2.14) 

The  superscripts  (1)  and  (2)  refer  to  the  half  plane  regions  to  the  left  and 
to  the  right  of  the  crack,  respectively.  The  singular  Integral  equation 
(2.13)  was  derived  in  [7]  using  integral  transform  techniques. 

The  governing  equation  for  this  problem  is  formulated  by  replacing  txxC 


In  eqn  (2.13)  by  the  difference  between  the  opening  pressure  and  the  spring 
induced  resisting  stress.  This  leads  to  the  following  equation: 


sin  c  j « 

3  (y* t) 

,  +  k(C0D)]  .  x.0,0<y<h  (2.15) 


As  stated  earlier,  the  opening  pressure  rxx  is  the  one  obtained  from  the  shear 
pullout  case,  and  k  is  the  equivalent  spring  constant  for  the  adhesive. 
Following  eqn  (2.14)  the  crack  opening  displacement  at  y  is  given  by: 


COO (y)  *  [u<2>  -  u^h  »  -/J  B(t)dt  ;  0  <  y  <  h  ,  x  «  0  (2.16) 

a  a  y  j 


Equations  (2.15)  and  (2.16)  are  solved  numerically  for  the  crack  opening 
displacements  and  the  dislocation  densities.  The  procedure  is  discussed  in 
detail  in  the  next  section. 


3.  Numerical  Analysis 


The  procedure  for  obtaining  numerical  results  for  the  shear  pullout  of 
the  rigid  insert  problem  is  discussed  first.  The  governing  integrodi fferen- 
tial  equation,  (2.6),  and  the  equation  of  equilibrium,  (2.10),  are  normalized 
by  introducing  the  following  variable  changes: 


y  a  ^  h(y  +  l)  ;  t  *  j  h(t  +  l) 


(3.1) 


The  function  D(t)  is  non-dimensional i zed  and  given  appropriate  square  root 
singularities  at  its  ends  by  making  the  following  substitution: 


0(t)  -  D(t)(l  -  l2)-0*5 


(3.2) 


With  these  changes,  equations  (2.6)  and  (2.10)  take  on  the  following  forms: 

.1  tHt)U--tV-5dt,  f  1  0(J)a.E2,-0.5  j(Jit-)d{ 
t  -  y 

,  1—  [0(y)(l-y2)~°,53  ;  x  -  0,  -1  <y,  1  <  1  ;  (3.3) 


>*sh  3y 


D(l)(l-l2)'0,5  dt  -  1  ;  x  *  0,  -1  <  l  < 


The  function  R(y,l)  is  as  given  next: 


(3.4) 


-l)2  «  2[2(£+l)-ic(y+E+2)]  a(y+l)(M) 


The  spring  constant  ks  Is  estimated  from  the  mechanical  properties  of  adhesive 


mate- idls  representative  of  the  ones  used  in  industry  [3].  A  simple  dimen¬ 
sional  analysis  carried  out  on  eqn  (2.1)  shows  that  the  appropriate  units  for 
the  equivalent  spring  constant  are  [_force/( length)^] .  If  G  is  the  shear 
modulus  of  the  adhesive  and  t  is  its  average  thickness  when  applied  to  such 
materials  as  aluminum  or  titanium  alloys,  then  k$  is  given  by: 

“s  *  |  (3.6) 

Numerical  results  were  obtained  for  several  values  of  ks.  The  range  of  values 
considered  was  based  on  two  adhesives  which  are  used  in  industry.  Their 
properties  will  be  discussed  in  the  next  section. 

The  kernel  of  the  right-hand  side  of  eqn  (3.3)  has  a  Cauchy-type  singu¬ 
larity.  The  numerical  method  introduced  by  Gerasoulis  and  Srivastav  [4],  [5] 
is  used  here  to  reduce  the  right-hand  side  of  (3.3)  to  a  system  of  linear 
algebraic  equations.  The  method  consists  of  replacing  the  inteyral  equation 
by  Integral  relations  at  a  set  of  points.  Piecewise  linear  functions  are  then 
used  to  approximate  the  Integrand  at  a  finite  set  of  (collocation)  points. 

The  values  of  the  unknown  function  at  those  points  are  then  obtained  via 
closed  form  Integration. 

Equation  (3.3)  Is  thus  reduced  to  a  system  of  linear  algebraic  equations 
which  may  be  written  In  matrix  form  as  follows: 

[A]  (0(f))  -  -  {*-  CD(y)(l-y2)’0-5])  (3.7) 

ci  ay 

Here,  the  symbol  [  ]  denotes  a  two-dimensional  matrix,  while  {  }  denotes  a 
column  matrix.  The  integration  points  ,  k*l,  2,  ...,  2N+1,  are  chosen  to  be 
equally  spaced  in  the  interval  -1  <  t,  _<  I,  and  the  collocation  points  y  are 
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2N.  The 


The  solution  is  obtained  using  an  iterative  method.  As  a  first  approxima¬ 
tion,  0( Ck )  are  set  equal  to  the  results  obtained  from  the  solution  of  the 
problem  which  assumes  that  the  insert  is  perfectly  bonded  to  the  half  plane, 
i.e.,  the  case  where  the  spring  stiffness  ks  tends  to  infinity.  Simple  matrix 
multiplication  leads  to  numerical  values  for  the  slopes  of  the  normalized 
stress  discontinuities  at  the  collocation  points  yj .  These  results  are  fitted 
using  cubic  splines,  and  then  numerically  Integrated  using  Gauss'  formula 
(see,  e.g.,  [10],  [11])  to  obtain  values  for  fl  at  These  values  are 
normalized  so  that  the  equilibrium  equation  (3.4)  is  satisfied,  and  the 
results  are  substituted  back  into  eqn  (3.7)  for  the  next  iteration.  This 
process  Is  repeated  until  a  convergence  criterion  Is  satisfied.  Upon  converg¬ 
ence,  the  results  are  used  to  compute  the  stresses  generated  in  the  contact 
region.  The  shear  stresses  are  obtained  from  eqn  (3.2).  To  compute  the 
normal  stresses,  equation  (2.11)  Is  first  normalized  using  eqns  (3.1)  and 
(3.2);  this  yields  the  following  expression; 


.  (ji  - p  <ru.  c/1 

xxx'y>  h  *(<+1)  'M  l_- 

+  0(£)(l-£2r°*5Il(y,f)dt}  ;  x  *  0,  -1  <  y  <  1 


(3.9) 


The  function  t(y,£)  is  given  by 


;  C(y.t) 


3<-l  +  2C3(t+l)-,c(y+l)3  8(y+l)(t+l)  j 
y+t*2  (y*t+2)2  (y+t*2)3 


(3.10) 


Equation  (3.9)  may  now  be  written  as  a  system  of  linear  algebraic  equations 
exactly  as  it  was  done  with  the  left  hand  side  of  eqn  (3.3).  Numerical  values 
for  the  normal  stresses  at  collocation  points  are  obtained  by  multiplying  the 
matrix  of  the  coefficients  of  this  system  of  equations  by  the  solution 
vector  (D(£)}. 

The  Gerasoul is-Sri vastav  technique  [4]  is  also  used  to  obtain  numerical 
results  for  the  opening  case.  First,  equation  (2.15)  is  normalized  and  the 
dislocation  densities  are  given  appropriate  singularities  at  the  ends  by 
making  the  following  substitutions: 


y  M  \  h(y+l)  ;  t  -  ^  h(t+l)  ;  (restated) 

S(t)  *  B(l)(l-E2r0,5 


With  these  substitutions,  equation  (2.15)  becomes: 


[  — ; —  *  ]  at  -  £  c-txx(o,J) 


(y*t+2)  (y+t+2)' 


+  k(^)]  ;  *  ■  0,  -1  <  y  <  1 


(3.1) 

(3.11) 


(3.12) 


The  equivalent  spring  constant  for  the  adhesive  In  tension,  k,  is  given  by  the 
following  relation: 


;v 

£ 

S: 


.N 

.*>■ 


i 


r.‘ 


7-’  V-  «■> 


k-|  (3.13) 

Numerical  results  were  obtained  for  values  of  E  corresponding  to  the  range  of 
G  values  used  in  the  shear  case.  A  detailed  discussion  of  this  is  presented 
in  the  results  and  discussion  section. 

Equation  (3.12)  is  reduced  to  a  system  of  linear  algebraic  equations  [4], 
and  then  placed  in  matrix  form  as  follows: 

CM]  Ml)}  =  c2{-rxx(0,y)  +  k(^)>  (3.14) 


The  matrix  [M]  contains  the  coefficients  of  the  unknown  normalized  dislocation 
densities,  and  c2  is  given  as 


*(<♦!)  h. 
2  P 


(3.15) 


Premultl plication  of  both  sides  of  eqn  (3.14)  by  the  inverse  of  [M]  leads  to 
the  following  relation: 


Ml)}  •  CgCM]-1  {-xxx(0,y)  +  k(^)> 


(3.16) 


The  solution  for  this  equation  Is  also  obtained  via  an  Iterative  method.  A 
first  estimate  of  the  values  In  the  column  matrix  {-Txx(0,y)  +  k(^-)>  Is  made 
by  setting  them  equal  to  some  percentage  of  the  opening  pressures  Txx(0,y). 

The  matrix  multiplication  Indicated  at  the  right-hand  side  of  eqn  (3.16)  leads 
to  the  Initial  estimate  for  the  dislocation  densities  B(tk).  These  values  are 
next  used  to  compute  the  corresponding  crack  opening  displacements  from  eqn 
(2.16).  First,  equation  (2.16)  is  normalized  by  making  the  substitutions 
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given  in  (3.1)  and  (3.11).  The  resulting  expression  may  be  written  in  the 
following  form: 


i  [COD]  *  \  [u[2)  -  u^]  »  -  t-  Jl  B(t)(l-t2)‘0,5dt  ;  -1  <  a  <  1 
c  a  c  x  Xa^u'a  - 

(3.17) 

The  values  ofl/2(C00)  are  needed  at  the  collocation  points  y  .  Values 
of  B(f)  have  been  calculated  at  the  integration  points  these  values  are 
numerically  interpolated  using  cubic  splines.  Then,  by  setting  a  *  yjt 
Gaussian  integration  leads  to  the  desired  values  of  V2  ( COD)  at  y  .  These 
results  are  then  substituted  into  the  right-hand  side  of  (3.16)  for  the  next 
iteration.  This  procedure  is  repeated  until  a  convergence  criterion  is 
satisfied. 


-14- 


a  "  » 


♦.  Results  and  Discussion 


The  numerical  analysis  was  carried  out  for  N  =  9  ana  N  =  17,  correspond¬ 
ing  to  19  and  35  Integration  points,  respectively.  There  was  no  significant 
difference  between  the  two  sets  of  results.  The  applied  load  P  and  the  depth 
h  were  set  equal  to  1  for  all  cases.  Results  were  obtained  for  a  range  of 
tensile  and  shear  elastic  moduli  E  and  G.  The  values  considered  were  set 
according  to  the  properties  of  two  commercially  available  adhesives.  The 
first  is  a  relatively  rigid  one,  known  as  FM-73,  manufactured  by  American 
Cyanlmld,  Bloomlngdale  Division.  It  is  an  elastomer  modified  epoxy 
material.  Its  properties  were  experimentally  determined  by  the  Kearfott 
Division  of  Singer  [12].  The  second  adhesive  material  is  of  relatively  low 
rigidity;  it  is  an  experimental  one  component  urethane  adhesive  produced  by 
Goodyear  coded  AX37J922  [13].  The  average  value  of  G  for  FM  73  is  given  in 
[12]  as  84,000  psl.  The  corresponding  value  for  the  Goodyear  urethane  system 
Is  15,000  psl  [13].  Average  values  for  the  installed  adhesive  thickness,  t, 
are  given  as  0.0047  Inches  (FM  73)  and  0.005  inches  (AX37J922).  Thus, 
following  eqn  (3.6)  the  corresponding  k$  values  are  1.789  x  107  lbs/1n3  and 
0.3  x  IQ7  lbs/1n3,  respectively. 

The  value  of  E  for  the  AX37J922  adhesive  Is  given  as  43,600  psl  [13]. 

This  makes  the  ratio  E/G  equal  to  2.907.  Poisson's  ratio  for  this  material  Is 
approximately  equal  to  0.45.  These  values  are  In  good  agreement  with  the 
theoretical  relation  between  E,  G,  and  v: 

E  •  2(l+v)G  (4.1) 


However,  the  experimentally  computed  value  of  E  for  FM  73  [12],  does  not 


correspond  to  the  value  given  for  G;  the  average  values  reported  are  E  3 
360,000  psi,  G  3  84,000  psi -  When  these  values  are  substituted  in  eqn  (4.1), 
the  resulting  value  for  v  is  1.14.  Since  v  cannot  be  greater  than  0.5,  the 
value  used  for  E  in  this  case  was  adjusted  to  be  consistent  with  the  one  for 
the  shear  modulus.  Taking  v  *  0.35  and  G  3  84,000  psi,  equation  (3.14)  gives 
E  3  226.00U  psi.  The  average  adhesive  thickness  t  is  given  as  0.0045  inches 
for  FM  73,  and  as  0.005  inches  for  AX37J922.  Thus,  the  corresponding  k  values 
are  5.022  x  107  lbs/in3  and  0.872  x  107  lbs/in3,  respectively.  Numerical 
results  were  obtained  for  four  Intermediate  sets  of  ks,  k  values  equally 
spaced  between  the  ones  already  given,  as  well  as  for  two  sets  of  k$,  k  values 
one  lower  than  the  AX37J922  ones,  and  the  other  higher  than  the  FM  73  ones. 

The  values  used  are  given  in  Table  1. 

Table  1.  Shear  and  Tensile  Spring  Constants 
Code  ks  (lbs/in3)  k  (lbs/in3) 

A  0.0022  x  107  0.042  x  107 

B  0.3000  x  107  0.872  x  107 

C  0.5878  x  IQ7  1.702  x  107 

0  0.8956  x  107  2.532  x  107 

E  1.1934  x  107  3.362  x  107 

F  1.4912  x  IQ7  4.192  x  107 

G  1.7890  x  107  5.022  x  107 

H  2.0868  x  107  5.852  x  107 

Two  adherend  materials  were  considered:  The  7075-T6  alumlnun  alloy  and  the 
T1-6A1-4V  titanium  alloy.  Both  of  these  materials  find  widespread  application 
in  aircraft  structures.  Material  constants  for  the  aluminum  alloy  are  taken 
as  E  •  10.4  x  10®  psi,  v  *  3.75  x  10®  psi,  and  v  =0.33;  the  ones  for  the 
titanium  are  £  3  16  x  10®  psi,  3  6 .4  x  10®  psi,  and  v  3  0.34.  All  results 


plotted  here  are  the  ones  obtained  for  the  plane  strain  case  with  N  =  17,  and 
with  T i -6A1-4V  being  the  adherend  material. 

The  normalized  shear  stress  distribution  generated  during  the  shear  pull¬ 
out  of  the  insert  is  plotted  in  Figure  2  for  all  ks  values  considered,  as  a 
function  of  the  (non-dimensionalized)  distance  from  the  free  surface.  As 
given  in  Table  1,  the  curve  Identified  by  A  corresponds  to  the  case  where  the 
adhesive  has  the  least  rigid  shear  spring  constant  k$  =  2.2  x  IQ4  lbs/in3.  At 
the  other  extreme,  curve  H  gives  the  shear  stress  distribution  for  the  case 
where  the  most  rigid  adhesive  is  used.  This  figure  shows  that  the  shear 
stresses  attain  their  maximum  values  near  the  tip  of  the  insert,  and  that  they 
are  smallest  near  the  free  surface  of  the  half  plane.  In  the  lower  limiting 
case.  A,  the  variation  is  not  detectable  and  the  stresses  are  constant  along 
the  insert.  The  difference  between  the  shear  stresses  near  the  tip  and  those 
near  the  free  surface  is  most  pronounced  for  the  case  involving  the  most  rigid 
adhesive.  As  the  rigidity  of  the  adhesive  increases  these  results  converge  to 
the  ones  obtained  for  the  “perfectly  bonded"  case  [6].  It  is  also  interesting 
to  note  that  the  curves  seem  to  "pivot"  about  a  point  near  which  the  value  of 
y  is  approximately  65X  of  h.  This  behavior  of  the  shear  stresses  is  predicted 
by  a  simple  1-0  model  analysis.* 

These  results  show  that  the  adhesive  will  tend  to  fall  (yield)  near  its 
tip  first.  This  will  result  in  a  redistribution  of  the  stresses  along  the 
bond  line.  This  type  of  non-linear  behavior  will  be  investigated  in  future 
work. 

The  non-dimensionalized  normal  stress  distribution  is  shown  in  Figure  3  as 
a  function  of  the  normalized  distance  from  the  free  surface.  For  clarity, 

*Th1s  1-0  model  comparison  was  pointed  out  by  Professor  P.J.  Torvik,  Air  Force 
Institute  of  Technology. 


only  the  two  extreme  cases  are  plotted.  These  figures  show  that  the  maximum 
normal  stresses  occur  near  the  free  surface,  and  they  tend  to  zero  as  y 
approaches  h.  It  is  also  observed  that  these  stresses  vary  significantly  with 
ks  only  at  their  extreme  values  near  y  =  0.  There,  the  stress  increases  with 
Increasing  ks.  Away  from  y  *  0,  the  normal  stress  distributions  seem  to  be 
independent  of  the  rigidity  of  the  adhesive. 

These  results  are  significant,  especially  when  viewed  in  light  of  the 
shear  stress  results  (Figure  2).  When  the  adhesive  fails  in  shear,  initially 
near  the  tip,  residual  compressive  stresses,  such  as  those  encountered  in 
fiber  reinforced  composite  materials  [1],  [14],  would  tend  to  keep  the 
debonded  surfaces  in  contact,  thus  inducing  friction  forces  in  the  failed 
region.  This  In  turn  would  tend  to  counter  shear  propagation  of  the  crack  at 
the  fiber-matrix  interface.  Since  the  normal  stresses  In  this  region  are 
negligible,  the  residual  stresses  would  dominate  and  the  crack  surfaces  would 
remain  In  contact.  However,  as  the  adhesive's  damaged  zone  increases  in 
length,  thus  moving  away  from  the  tip,  the  normal  stresses  increase,  and  at 
some  distance  from  the  tip  they  will  balance  the  residual  compressive 
stresses,  thus  eliminating  friction.  Beyond  that  point,  the  growth  of  the 
damaged  zone  may  accelerate  under  the  combined  action  of  the  shear  stresses 
and  of  the  opening  pressure,  until  complete  debonding  is  reached.  This  sug¬ 
gested  failure  process  seems  to  be  In  very  good  agreement  with  the  experi¬ 
mental  results  obtained  by  Atkinson,  et  al.  [15]  for  the  axlsymmetric  case. 

The  results  obtained  from  the  solution  of  the  problem  involving  the  edge 
crack  subjected  to  opening  pressure  are  plotted  in  Figure  4.  The  crack  open¬ 
ing  displacements  are  symmetric  about  the  x  *  0  axis  (see  Fig  1).  The  curves 
shown  represent  the  opening  at  either  side  of  the  insert  as  a  function  of 
(normalized)  distance  from  the  free  surface.  Thus,  these  are  the  displace- 


ments  which  are  multiplied  by  the  tension  spring  constant,  k,  to  obtain  the 
resistance  to  opening  due  to  the  adhesive. 

These  results  show  that  the  opening  displacements  increase  with  decreasing 
adhesive  stiffness,  k.  This  is  what  one  would  expect,  considering  the  opening 
pressure,  txx,  does  not  vary  substantially  for  different  adhesive  materials 
(ks),  as  shown  in  Figure  3,  while  the  resistance  to  opening,  k( COO/2) , 
increases  considerably  with  increasing  k.  Thus  the  magnitude  of  the  right- 
hand  side  of  eqn  (2.15)  is  very  sensitive  to  the  magnitude  of  k.  This  is 
shown  in  Figure  5,  where  the  "net"  opening  pressure,  -[-txx  +  k(-^p)],  is 
plotted  for  all  cases  as  a  function  of  distance  from  the  free  surface.  This 
figure  shows  good  agreement  between  the  applied  “net"  opening  pressures 
(increasing  with  decreasing  k)  and  the  crack  opening  displacements  shown  in 
Figure  4. 

It  should  also  be  noted  that  the  crack  opening  displacements  near  the  tip 
(y  ■  h)  are  between  one  and  two  orders  of  magnitude  smaller  than  the  ones  near 
the  surface  (y  *  0).  This  behavior  Is  consistent  with  the  variation  in  the 
net  opening  stresses  plotted  In  Figure  5.  Thus,  the  crack  is  essentially 
closed  over  a  portion  of  Its  length  near  the  tip  of  the  insert,  and  the  length 
of  the  closed  segment  Increases  with  k.  It  should  be  emphasized  here  that 
this  analysis  does  not  account  for  any  residual  compressive  stresses  that  may 
be  Introduced  during  the  manufacturing  processes.  The  presence  of  any  such 
stresses  would  clearly  tend  to  Increase  the  length  over  which  the  crack  Is 
closed. 
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